Superconductivity in the Kondo lattice: a mean-field approach 
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We calculate the superconducting critical temperature Tc and the Kondo temperature Tk of 
the Kondo lattice, decoupling the Kondo exchange interaction in the unrestricted Hartree-Fock 
(HF) Bardeen-Cooper-Schrieffer (BCS) approximation. We obtain that both Tk and Tc have an 
exponential dependence in the Kondo coupling Jk- For optimum doping and realistic parameters, 
both temperatures fall in the experimentally observed range. 
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I. INTRODUCTION 

One |©f the main characteristics of heavy-fermion 
systemstl is the competitiou j-between Kondo effect and 
magnetic long-range order,B'0 both stemming from the 
local antiferromagnetic (Kondo) coupling between lo- 
calized moments and conduction-electron spins. The 
non-magnetic compounds, like CeAls, present a low- 
temperature state resembling a Fermi liquid with en- 
hanced carrier mass, the so called heavy-fermion regime. 
In addition to that, supercondtuctivity is observed be- 
low Tc ~ IK in some systems,Qi3 the best known being 
CeCu2Si2, UBei3, and UPta. Although the nature of su- 
perconductivity in these systems is still unclear, it has 
been suggested that the Kondo interaction itself might 
provide a pairing mechanism.cl Indeed, there is experi- 
mental evidence that heavy electrons are involved in the 
superconducting state, as a high value is observed for the 
specific-heat jump at the superconducting transition. 

Lately, the competition between Kondo effect and 
magnetism has been studied in a mean-field approxi- 
mation based xin a simple decoupling of the Kondo ex- 
change termJj'El but superconductivity has not been stud- 
ied within a similar approach. If fluctuations are not 
included, this mean-field approximation is equivalent to 
considering only part of the terms of an unrestricted 
Hartree-Fock (HF) decoupling of the exchange interac- 
tion. In this work, we include all relevant terms in this 
HF decoupling. Among them, a singlet pairing coupling 
appears, which we treat in the Bardeen-Cooper-Schricffer 
(BCS) approximation, leading to superconductivity in 
the model. Although there are experimental grounds to 
the notion that superconductivity in heavy-fermion sys- 
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tems presents non-trivial symmetry properties,u and can 
even coexist with magnetic ordering, it is instructive to 
explore, in a simple theoretical approach, the existence 
of a superconducting solution in which the pairing mech- 
anism is the same local Kondo interaction that is also 
responsible for magnetic ordering and compensation of 
local moments. 



II. MODEL AND MEAN-FIELD 
APPROXIMATION 



We consider the standard Kondo-lattice model used to 
describe heavy-fermion systems, writing the Hamiltonian 
in the form 



H = -t 



{cl^Cja + H.c.) + Jk ^ Si • Si 



(1) 



creates a conduction electron with spin a at site 
(si) is the localized (conduction) spin operator 



where cl^ 
i, and Si 
at site i. 

The|-fIF-BCS approximation has been described 
before)13 and has been used, for example, to study d-wa^ 
superconductivity in a generalized Hubbard modcl,Ell 
with results which are in qualitative-agreement with ex- 
act diagonalization of finite systems.E3 We use a fermionic 
representation of the localized spins, assigning them to / 
electrons, whose occupation number is restricted to one 
per site. Then, neglecting expectation values which do 
not conserve spin (e.g, triplet pairing), the unrestricted 
HF decoupling of the interaction term reads 



1 



l)(4 



H.c. 



2 (/it/i 



(2) 
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where a = —a. The constraint of one locaHzed (/) par- 
ticle per site is taken only on average, i.e., 



(3) 



N being the number of sites, and is imposed by adding 
to the Hamiltonian the term 



f 2-^ ficrfic^' 



(4) 



where Ef is a Lagrange multiplier. This procedure re- 
duces the Hamiltonian to an effective one-body problem. 

In the following we consider only the non-magnetic ho- 
mogeneous solution with singlet superconductivity. The 
order parameters are A — (cl^fia), characterizing the sin- 
glet Kondo state,El and t] ~ {fi^^l^), which is the usual 
superconducting order parameter. Both can be made real 
through a gauge transformation. 

In absence of superconductivity (rj = 0), the problem 
is easily solved. The electronic structure oi H + He is 
composed of two hybrid bands with energies 



Ef^ = 'J^±^ - (+ 




(5) 



where is the dispersion relation of the unperturbed 
conduction band, and 



V 



(6) 



Both, Ef and A must be determined selfconsistently. The 
latter shows the usual behavior of a mean-field order pa- 
rameter as a function of temperature, going to zero at 
a critical temperature which we identify as the (mean- 
field) Kondo temperature Tk. In the limit A ^ 0"*", the 
self consistent conditions lead to a simple equation for 
Tk- 



3 tanh(S^) 



8N 



(7) 



The chemical potential /i in this limit is determined by 
the number of conduction electrons of the unperturbed 
band. 



(8) 



with /(e) standing for the Fermi function, and Pq^e) be- 
ing the unperturbed density of states (DOS) per spin. In 
the following we take, for simplicity, a parabolic DOS, 
Poie) — 3(1 — e^)/4, choosing the half band width of the 
unperturbed band as the unit of energy. 

Since H is electron-hole symmetric, we can assume, 
without loss of generality, Uc < 1, so that the Fermi 
level falls in the lower hybrid band. The pairing terms 
in H can be written in terms of the hybrid operators 



describing the quasi-particles in the lower band, and we 
can safely neglect the upper one, as the superconducting 
critical temperature Tc turns out to be much lower than 
Tk. After some algebra, we arrive at the following four 
equations, which determine Ef, A, p, and the supercon- 
ducting critical temperature Tc'. 



Hf 



K 



^/+"= = |E[/(^K) + /(i^K)], 
K 



1 = 



1 



167V 2^rl{E^-p) 



tanh ' ^ 



Tr_ 



(9) 



where = 1/2 + {Ef - eK)V(4rK), Yl = l- X^, and 
rK = [iEf-eK)V4 + V^]^/\ 

III. RESULTS 

In Fig. 1 we present Tk and the superconducting criti- 
cal temperature Tc as a function of conduction-band fill- 
ing Tic for a typical value of Jk . The Kondo temperature 
Tk increases with tIc and has its maximum at ric = 1 
(as explained above, TK(nc) = Tk(2 — ric)). Instead, Tc 
should vanish at ric = 1 because this case corresponds 
to a Kondo insulator, with the Fermi level falling in a 
gap whose magnitude is larger than the effective pairing 
interaction. In other words, while the effective pairing 
interaction and the density of states increase as Uc tend 
to one, the effective number of carriers around the Fermi 
level tends to zero. As a consequence of this competition, 
Tc shows a maximum around ric 0.65-0.7, whose posi- 
tion is not strongly dependent on the coupling, as shown 
in Fig. 2. Assuming that the magnitude of the unper- 
turbed half band width is of the order of 1 eV, then, for 
Jk ~ 0.3 eV we obtain Tc ~ 1 K at optimum doping (see 
Fig. 1). 




FIG. 1. Kondo temperature Tk and superconducting criti- 
cal temperature Tc as a function of doping for Jk = 0.3. 
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FIG. 2. Superconducting critical temperature Tc as a func- 
tion of band filling for two different values of the Kondo cou- 
pling. 
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FIG. 3. Kondo temperature Tk and superconducting criti- 
cal temperature Tc as functions of Jk for ric — 0.7. 



Figure 3 shows Tk and Tc as functions of Jk for 
ric = 0.7 (near optimum doping). The dependence of 
both temperatures with Jk can be fitted very well by 
exponential functions: 



Tk = Acxp(-a/JK), 
rc = Sexp(-/3/JK), 



(10) 



as clearly evidenced by the logarithmic plots of Fig. 4. 
The meaning of the coefficients is not obvious. The pa- 
rameter values used in Figs. 2 and 3 yield A — 0.7289, 
a = 1.876, B = 0.3250, and /? = 2.425. For a single 
Kondo impurity, in the limit of infinite band, width, it is 
known that T™'' = (1/po) 6xp[— l/(po JK)],li3 where Pq is 
the magnitude of the unperturbed density of states, as- 
sumed constant. The natural extension to the case of a 
finite band width would be to write the single-impurity 
Kondo temperature as the first of equations (M) with 



a™p = 1/Po(sf)' which would have the value 1.3905 for 
the example considered. However, in the present case 
the coefiicients are not only affected by the finite band 
width but also by the mean-field approximation and the 
fact that we are dealing with a Kondo lattice instead of 
a single impurity. The dependence of Tc with Jk is even 
more difficult to predict. From usual BCS theory one 
might expect an exponential dependence on the product 
of the interacting density of states and the effective pair- 
ing interaction, if the latter is small. This in turn, is 
proportional to Jk and (see Eq.(||)), with a prefactor 
that, being generally of the order of magnitude of the 
band width, should vanish as nc 1. 



IV. DISCUSSION 

We have shown that the Kondo-lattice model, treated 
within the HF BCS approximation, can lead to a su- 
perconducting phase within the Kondo (heavy-fermion) 
regime, with the Kondo temperature Tk and supercon- 
ducting critical temperature Tc of the order of magni- 
tude observed in experiments. While both temperatures 



depend exponentially on the exchange interaction, their 
doping dependences are different, and Tc vanishes for a 
half-filled band. If one wants to make quantitative com- 
parisons with experimental results, an appropriate choice 
of parameters that gives a good agreement for Tc will 
yield a slightly overestimated Tk. This is typical of a 
mean-field solution, and one expects the true Kondo tem- 
perature to be reduced by fiuctuation effects. Actually, 
we would like to remark that Tk signals a second-order 
phase transition in the mean-field approach, while only 
broad features characteristic of a crossover behavior are 
observed in experiments. 

It is interesting to notice that other approximations 
can be obtained if one explicitly uses the constraint 
ni = 1 to rewrite parts of the Hamiltonian before per- 
forming the HF decoupling. For example nfn? = n1 
ceases to be an exact equality once the first term is de- 
coupled. The use of different approximations affect the 
strength of the different decoupled terms, and therefore 
affect the relative stability of the different phases. How- 
ever, assuming that magnetic ordering does not occur, 
the qualitative behavior of Tk and Tc should be the same 
as in our approach. 
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FIG. 4. Logarithmic plots of Tk and Tc as functions of 
1/ Jk for ric = 0.7, showing the straight lines that fit Eqs. (^o|) . 
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We have not analyzed here the stabihty of our solution 
with Tjespect to magnetic ordering. However, previous 
resultaa have shown that within a mean-field approach 
one needs an explicit intersite exchange in order to ob- 
tain magnetic long-range ordering. Even in that case, the 
non-magnetic solution should be stable for moderate or 
large Kondo interaction Jk. On the other hand, the pos- 
sibility of obtaining a wave-vector dependent supescon- 
ducting gap, as suggested by experimental resultsja has 
more to do with the model itself, as one should consider a 
non-local Kondo interaction. The present approach can 
be easily extended to such a case. 
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